Electron dynamics in crystalline semiconductors 
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Electron dynamics in crystalline semiconductors is described by distinguishing between an in- 
stantaneous velocity related to electron's momentum and an average velocity related to its quasi- 
momentum in a periodic potential. It is shown that the electron velocity used in the theory of 
electron transport and free-carrier optics is the average electron velocity, not the instantaneous ve- 
locity. An effective mass of charge carriers in solids is considered and it is demonstrated that, in 
contrast to the "acceleration" mass introduced in textbooks, it is a "velocity" mass relating car- 
rier velocity to its quasi-momentum that is a much more useful physical quantity. Among other 
advantages, the velocity mass is a scalar for spherical but nonparabolic energy bands e(fc), whereas 
the acceleration mass is not a scalar. Important applications of the velocity mass are indicated. A 
two-band k ■ p model is introduced as the simplest example of a band structure that still keeps track 
of the periodic lattice potential. It is remarked that the two-band model, adequately describing 
narrow-gap semiconductors (including zero-gap graphene), strongly resembles the special theory of 
relativity. Instructive examples of the "semi-relativistic" analogy are given. The presentation has 
both scientific and pedagogical aspects. 

PACS numbers: 71.28.+d, 73.61. Ey, 72.20.-i 



I. INTRODUCTION 



The phenomenon of Zitterbewegung (ZB, trembling 
motion), devised in 1930 by Erwin Schrodinger [l|, has 
been for the last 80 years a subject of controversy and ex- 
citement. The interest in this phenomenon experienced 
a strong revival in 2005, when it was demonstrated that 
the trembling motion can occur also in solids @, Q . Since 
then there has been a real surge of papers proposing ZB in 
various periodic systems, as reviewed in jjj. The nature 
of ZB in solids was investigated and it was shown that, 
in its "classical" form analogous to ZB in a vacuum [![, 
the trembling motion represents oscillations of velocity 
when an electron moves in a periodic potential of the lat- 
tice The situation resembles a roller-coaster: when 
the train moves upwards gaining the potential energy, it 
slows down; when the train goes down losing the poten- 
tial energy, it accelerates. However, in the solid state 
literature the largely prevailing picture is based on the 
Bloch theorem in which electrons are treated as quasi-free 
particles with a modified (effective) mass. This picture 
suggests that the electrons move in a solid with a constant 
velocity. The same approach is used in the transport 
theory, in which carrier velocity is assumed constant and 
equal to v — hk /m* , where k is the wave vector and m* 
is the effective mass. To the author's knowledge, there 
exist only two textbooks discussing an instantaneous car- 
rier velocity in a crystal: "Wave Mechanics of Crystalline 
Solids" by R. A. Smith @ and "Semiconductor Physics" 
by P. S. Kireev Q. 

The problem arises: how to reconcile the two pictures? 
This is the first purpose of our work. We show that 



the above question is related to the difference between 
carrier's momentum and quasi-momentum in a periodic 
potential. Our second purpose is to introduce properly 
an effective mass of carriers. In solid state textbooks 
the effective mass is always defined as a quantity relat- 
ing an external force to carrier's acceleration. We show 
that it is by far more useful to define the effective mass 
as a quantity relating the average carrier velocity to the 
quasi-momentum hk. The "velocity mass" is a scalar for 
spherical nonparabolic energy bands e(k), whereas the 
"acceleration mass" is not. Important applications of 
the velocity mass are indicated. In addition, we briefly 
describe a "semi-relativistic" behavior of charge carri- 
ers in narrow-gap semiconductors including monolayer 
graphene. This feature was discussed in the past but 
gained new significance with "the rise of graphene" Q 
and the advancement of Zitterbewegung. To ensure the 
completeness and continuity of presentation, we include 
in our exposition a few elements which are already known 
from the literature. As it stands, the text has both sci- 
entific and pedagogical aspects. 



II. ELECTRONS IN A PERIODIC POTENTIAL 

We begin by general considerations concerned with the 
motion of charge carriers in crystalline solids. The Hamil- 
tonian for an electron in a periodic potential V(r) is 



H 



P 

2ni{) 



V(r), 



(1) 



where m is the free electron mass. The periodicity sig- 
nifies V(r) — V(r + a) for a being a lattice vector. The 
velocity operator is given by the Hamilton equation 
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The same result is obtained from the relation ii 
dxi/dt — (l/iti)[Xi, H]. The acceleration operator is 



at ihm 



--^ = -ir, (3) 

mg oxi mo 



where Ff r = —dV/d£i is a periodic force acting on the 
electron moving in a periodic potential. Equation (3) is 
equivalent to the second Newton law of motion in an op- 
erator form. It follows from Eq. (3) that the momentum 
operator does not commute with the Hamiltonian (1), so 
it is not a constant of the motion. In consequence, the 
velocity of Eq. (2) and the acceleration of Eq. (3) are 
also not constants of the motion. It is intuitively clear 
that, since the potential and the resulting force in Eq. (3) 
are periodic in r, the acceleration and the velocity will 
also be periodic functions of r. This result has an el- 
ementary classical interpretation. Classically, the total 



electron energy is e 



72 + V{r), so if the potential 



energy oscillates, the kinetic energy (i.e., the velocity) 
also oscillates to keep the total energy constant. This 
is, in fact, the physical origin of the trembling motion of 
electrons in crystalline solids, see Ref. The results 
given in Eqs. (2) and (3) apply to any Hamiltonian with 
a scalar potential. The specificity of a crystalline solid 
is that the potential is periodic, so the Bloch theorem 
applies. Thus 



Hij) nk {r) = e n (k)ip nk (r), 



(4) 



where e n (k) is the energy of the n-th band depending 
on the wave vector k. The Bloch state is ipnk{r) = 
exp(ifc ■ r)u nk (r) in which the Bloch amplitude has the 
same periodicity as the potential in Eq. (1), i.e. u nk (r) = 
u n k(f + o)- The quantity hk is an eigenvalue of the 
quasi-momentum operator P, which should be carefully 
distinguished from the standard momentum operator p 
introduced in Eq. (1). The Hamiltonian (1) is invariant 
with respect to transformations possessing the symme- 
try of the potential and there should exist a constant of 
the motion corresponding to this invariance. This con- 
stant of the motion is precisely hk. It means that the 
Bloch state ip n k{r) should also be an eigenstate of the 
quasi-momentum P, i.e. there should be by definition 



Pipnk{r) = hktp nk (r) 



(5) 



We try to find an explicit expression for P looking for 
the quasi-momentum operator in the form (see Ref. |7|) 

P = p + ihj(r), (6) 

in which -f(r) is a function of coordinates. We have 

Ptpnk = hktp nk + ih-y(r)ijj nk - ihe lkr V r u nk 

= hkip n k + ih[j(r) - V r (\nu nk )]ip nk . (7) 

Putting 7(r) = V r (lnu„ fc ) we get Pip nk = hkip nk , so 
that Eq. (5) is satisfied if 



In the second term in Eqs. (7) and (8) the differentiation 
acts only on the expression in parentheses. Equation (8) 
is instructive, as it shows explicitly that the operators of 
momentum and quasi-momentum are distinctly different. 
By using Eqs. (4) and (5) one easily shows that HP = 
PH, so that the quasi-momentum hk is really a constant 
of the motion. To say it differently 



(9) 



which means that the periodic potential V(r) in Eq. (1) 
does not change the quasi-momentum P whereas, as fol- 
lows from Eq. (3), it periodically changes the momen- 
tum p and velocity v. Still, the electron does not radiate 
because it is in the Bloch eigenenergy state. Suppose 
now that, in addition to the periodic potential V(r), 
the electron experiences an additional nonperiodic po- 
tential U ex (r). This potential can be due to an external 
field, an impurity, a defect, etc. Then the total Hamilto- 
nian is 



P 

2m 



+ V{r) + U ex (r). 



It is easy to see that 

dpi 
dt 

dt 



= Ff 



(10) 

(11) 
(12) 



where F ex = —~V r U ex (r). Thus the momentum is 
changed by both periodic and nonperiodic potentials, 
whereas the quasi-momentum is changed only by the ad- 
ditional nonperiodic potential. 

The question arises how to reconcile the oscillating 
electron velocity v(t) described in Eqs. (2) and (3) with 
the velocity appearing, for example, in the transport phe- 
nomena and other kinetic effects in crystalline solids. The 
time-dependent instantaneous velocity is given in general 
in the Heisenberg picture by 



v(t) — exp(iHt/h)v exp(—iHt/h), 



(13) 



and it is this velocity operator that leads to the trembling 
motion, see Ref. Let us calculate an average of v(t) 
on the Bloch state ipnk(r). We obtain by a simple ma- 
nipulation 

v = (^ak\e lflt/h ve-^ h \^ nk ) 

= (^nk\e te ^ t/h ve- ie ^ h \i; nk ) = (ip nk \v\^ nh ),(U) 

where v = p/rrio. Thus the average velocity in the Bloch 
state is time independent because of the basic property 
of Eq. (4) . The average velocity has been calculated in 
various ways, see Refs. [9l4ll|. Below we use the method 
based on the Hellmann-Feynman theorem [T3|. Let us 
first write the Schrodinger equation for the Bloch ampli- 
tude. As follows from Eq. (4) 



-ihV r + ih"V r (hiu nk ) 



(8) 



H u (r,p;k)u nk (r) = e n (k)u nk (r), 



(15) 
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where 



H u (r,p;k) = — (p + hk) 2 + V(r), (16) 
depends parametrically on k. We have 



V fc e n (fe) = {u nk \V k H u \u nk ) = (u nk \ (p + hk)\u nk ) 

mo 

h 2 

= — (u n k\ - i~V r + k\u nk ), (17) 
m 

in which the first equality follows from the Hcllmann- 
Feynman theorem. Further 



(-iV r + k)u nk = (-iV r + k)e ik ' r tp nk 



(18) 



so that 



-V fe e„(fe) = -i — (u nk \e ik ' r V r \i> nk ) 
h too 

P 

= {lpnk\ \tpnk) = V n (k). (19) 

m 

The result (19) is simple and important. It relates 
the average of instantaneous velocity in the Bloch 
state {ipnk\p/ rr i'o\' l l J nk) to the electron energy e n (k) given 
as a function of the quasi- momentum hk. Below we will 
consider a specific energy band, so we drop the band in- 
dex n. 

Now we want to associate the above results with an 
effective mass of charge carriers in an energy band. In 
textbooks one considers standard parabolic and spherical 
energy bands described by the energy-wave vector rela- 
tion e = ft 2 fc 2 /2mQ, where is a constant effective mass 
at the band edge. It is then shown that such a mass re- 
lates carrier's acceleration to an external force. However, 
we want to consider a more general case of spherical but 
nonparabolic energy bands in which the energy depends 
on the absolute value of the wave vector in an arbitrary 
way, i.e. e = e(k). In fact, many III-V semiconducting 
compounds (InSb, InAs, GaSb, GaAs, InP) as well as II- 
VI compounds (HgTe, CdTe, HgCdTe, HgSe) and their 
alloys possess the conduction bands of this type. In con- 
trast to the procedure adopted in textbooks, we define an 
effective mass not by a relation between an external force 
and acceleration, but as a quantity relating the average 
velocity v to the quasi- momentum hk. Thus we define 
the effective mass by the equality 



m*v = hk, 



(20) 



where v is given by Eq. (14). Since v and hk are vectors, 
the mass to* is in principle a 3 x 3 tensor. Using Eq. (19) 
and the sphericity of the band we calculate 



de de dk de k\ de 1 
hdh hdk dki hdk k hdk k 



where in the last term we adopt the sum convention over 
the repeated coordinate subscript j = 1,2,3. Using the 
definition (20), the inverse mass tensor is 



r, = | — ] hkj. 



By equating Eq. (21) with Eq. (22) we obtain 



1 

m* 



de 1 



h 2 dkk 13 ' 



(22) 



(23) 



Thus the inverse mass tensor is a scalar for a spherical 
energy band 



1 1 de 

to* h 2 k dk 

The average velocity is finally, see Eq. (22), 

hk 

v = . 

to* 

Recalling that v = p/mo, see Eq. (19), we can write 



(24) 



(25) 



P 

to 



hk 

TO* 



(26) 



Sijkj, (21) 



which shows an analogy between the average momentum 
in the Bloch state tp k and the quasi-momentum. How- 
ever, it is known that for electrons and light holes in semi- 
conductors there is usually to S> to*, so that p ^> hk. 
This shows once again the difference between momentum 
and quasi-momentum. 

Equation (25) represents the basic formula for velocity 
used in the description of charge carriers in semiconduc- 
tors and metals. Here we have obtained this formula 
with two important qualifications. First, on the left- 
hand side we have the average velocity of a carrier in 
the Bloch state, not the instantaneous velocity consid- 
ered in the beginning, see Eq. (2). Second, on the right- 
hand side we have the velocity effective mass defined in 
Eq. (20). This mass depends in general on carrier's en- 
ergy (or wave vector). Since the average velocity v is 
expressed by the first derivative de/dk, the velocity effec- 
tive mass to* is also related to the first derivative de/dk. 
On the other hand, the " acceleration" effective mass My 
relating force to acceleration, as introduced in textbooks, 
is given by the second derivative of energy with respect 
to fe, so this mass does not enter into the basic for- 
mula (25), unless one takes the simplest energy band 
described by e = h 2 k 2 /2rriQ. As is easy to see, in this 
particular case both masses are equal to TOg . We empha- 
size that the velocity mass, defined in Eq. (20), is much 
more useful than the acceleration mass defined in text- 
books. In particular, it is the velocity mass that defines 
carrier's mobility and is measured in various experiments. 
We discuss this point below. 

To conclude this section, we relate the first deriva- 
tive de/dki to the group velocity of a carrier in a periodic 



4 



potential. Let us form a wave packet of Bloch states 
f(r) = J a(ky kr u k (r)d 3 k 

~ u ko (r) [ a{k)e ik - r d 3 k. (27) 



It is assumed that the packet is narrow in k space and 
it is centered around the value of fco. This means that 
in the coordinate space the packet extends over several 
unit cells. Thus the amplitudes a(k) in Eq. (27) are non- 
vanishing only for small values of q = k — fc . In conse- 
quence, it is possible to take an average value of u ko (r) 
out of the integral sign. We further have 



f(r) = u ko (r)e 



iko-r 



a(ky qr d 3 q. 



(28) 



The integrand on the right-hand side of Eq. (28) is for- 
mally identical to that of a free particle and one may 
apply to this wave packet the well known arguments de- 
termining the group velocity, which gives 



duj 



hdh 



(29) 



where the derivative is taken at k = fco. Thus the average 
velocity given in Eq. (19) is also the group velocity of the 
carrier. 



III. TWO-BAND MODEL. SEMIRELATIVITY. 
GRAPHENE. 

Now we consider an instructive and sufficiently gen- 
eral example of a band structure in semiconductors in 
order to illustrate consequences of the above formalism. 
As mentioned in relation to Eqs. (1) and (4), the Bloch 
states are solutions to the eigenenergy equation with the 
Hamiltonian having a periodic potential of the crystal 
lattice. However, it is known that, for questions related 
to the band structure near a specific point of the Bril- 
louin zone or to problems of carriers in external fields, it 
is more practical to work with the Luttinger-Kohn (LK) 
representation (see Refs. [l3| and [IH). The LK func- 



tions are Xnk{r) 



' r u n o(r), where u n o(r) are the 



Bloch periodic amplitudes taken at a fixed point fco of 
the Brillouin zone. We take for simplicity fco = 0, i.e. 
the zone center. It is clear that the LK amplitudes sat- 
isfy the eigenenergy equation 



V 

2mo 



V(r) 



U n = EnOWnO, 



(30) 



where e n0 is the energy of the n-th band at fc = 0. One 
can show that the LK functions form a complete orthog- 
onal set, so one can represent a Bloch state as 



i'feW = ^^c:'(fc)u„ (r) 



(31) 



in which c™ (fc) are fc-dependent coefficients. The sum- 
mation is over all bands n. What follows is the stan- 
dard procedure of transforming a differential eigenvalue 
equation into an algebraic problem. By inserting the 
form (31) into initial Eq. (1), using Eq. (30), multiplying 
on the left by u n 'o and integrating over the unit cell, one 
obtains 



E 



{-(' + e n o)S n >„ H fc • p n i n 

m 



0, 



(32) 



for n' = 1,2,3.... Here e' = e — H 2 k 2 /2mo and p n 'n = 
(u n o\p\u n o) are the interband matrix elements of momen- 
tum. Equation (32) represents an infinite set of equations 
for c™ coefficients and the condition of non-trivial solu- 
tions determines the energies e n >(k). We now assume 
that an energy gap e g between the conduction and va- 
lence bands is much smaller than other gaps of inter- 
est, so we can neglect the distant bands and keep in 
Eq. (32) only the two close bands. In addition, we ne- 
glect the free electron term h 2 k 2 /2mo in the energy as 
it is small compared to the effective mass term, see be- 
low. Taking the zero of energy in the middle of the gap, 
so that eio = +e 5 /2 and e 2 o = — e g/2, the set (32) is 
reduced to 



+e g /2 7Ti2 • hk 
7T 2 1 • hk —e g /2 



(33) 



where 7Ti2 = PyiItciq and similarly for 7T2i. Solving the 
above set for the energies one obtains 



e(k) = ± 



h 2 k 2 

' 2m* 



1/2 



(34) 



if we assume the simplest symmetry of the matrix ele- 
ments giving 27Ti27ri2/e 5 = (l/m^Sij. Here is the 
electron effective mass at the band edge. Plus and minus 
signs correspond to the conduction and valence bands, 
respectively. Bands described by Eq. (34) are spherical 
and nonparabolic. For H 2 k 2 /2rriQ <C e g /2 one can expand 
the square root and obtain e(k) = ±(e g /2 + h 2 k 2 /2mg), 
so that for small k values the bands are parabolic, while 
for large k values they are linear in fc. Using for e(fc) 
relation (34) one can easily calculate the energy depen- 
dence of the velocity mass m* given by Eq. (24). For the 
conduction band one obtains 



2e 



m 



nin 



(35) 



'0; 



as it 



At the band edge e = +e g /2 there is m* 
should be. The band-edge mass ttiq in most semicon- 
ducting materials is much smaller than the free electron 
mass mo, so neglecting the free electron term h 2 k 2 /2mo 
in Eq. (32) was justified. 

It was remarked that the two-band fc • p model (2BM) 
for the band structure of semiconductors closely resem- 
bles the description of free relativistic electrons in a vac- 
uum p^| - [l7| . The Hamiltonian (33), having the quasi- 
momentum terms off the diagonal, looks very much 



like the Dirac equation without spin, while the disper- 
sion (34) is analogous to the relativistic relation E = 
±y/ (moc 2 ) 2 + c 2 p 2 with the correspondence p — > hk and 



2moc — > €„ 



m 



(36) 



It is easy to determine the maximum velocity u in the 
2BM 



2moc 2 
2ma 



1/2 



2m? 



1/2 



(37) 



In light of our previous considerations, u is the maximum 
average velocity in the Bloch state. The value of u can 
be determined by measuring the energy gap e g and the 
band-edge mass in a semiconductor material. It turns 
out that the velocity u is almost the same in different ma- 
terials and is given by u « 10 8 cm/s, i.e. it is about 300 
times smaller than the maximum velocity for relativistic 
electrons in a vacuum c. Using Eq. (37) for u one can 
rewrite Eq. (35) in the form 

e = m*u 2 . (38) 

This is equivalent to the famous Einstein formula: E — 
mc 2 relating the energy to the mass. The Compton wave- 
length A c = h/rriQC, playing an important role in the rel- 
ativistic quantum mechanics, also has a corresponding 
length in the two-band k ■ p model, see Ref. [2] 



ft 



= ft 



1/2 



(39) 



This length determines the amplitude of Zitterbewegung 
oscillations mentioned in the Introduction, see In 
narrow-gap semiconductors one can have mj = 5 x 
10~ 2 77i and, since c s» 300u, one obtains Xz ~ 2 x 
10 4 A C « 50A, i.e. a sizable length for nanostructures. 
The dispersion relation (34) can be rewritten in terms 
of u and Xz in the form 



;(fc) = ±hu (X z 2 + k 2 ) 



1/2 



(40) 



For k 2 > 0, Eq. (40) describes the conduction and 
light-hole bands. For k 2 < 0, that is for imaginary val- 
ues of k — in, this equation describes the dispersion 
in the gap. The latter can be determined in metal- 
semiconductor tunnelling experiments. In Fig. 1 we show 
the results of Parker and Mead [l8[ for InAs, as de- 
scribed by Eq. (40) with the use of adjustable parame- 
ters u and Xz ■ One obtains a very good description which 
confirms the validity of the two-band model for narrow- 
gap materials. In particular, one obtains the above men- 
tioned large value of Xz ■ 

Regarding the phenomenon of Zitterbewegung we want 
to emphasize the following subtle point. If one calcu- 
lates the velocity operator using the matrix Hamilto- 
nian (33): f), = dH/dhki, the velocity matrix does not 
commute with the Hamiltonian: iiiH — Hii ^ 0, so the 




0.2 



- 0.1 



0.0 §- 

LU 



-0.1 



0.2 



5 4 3 2 
-k 2 x 10 4 (A" 2 ) 



FIG. 1: Energy- wave vector dependence in the forbidden 
gap of InAs. Various symbols show experimental data of 
Parker and Mead ["HI ], the solid line is theoretical fit using 
Eq. (40). The determined parameters are Xz = 41.5 A and u 
= 1.33x10 s cm/s. After Ref. 0. 



velocity depends on time and it is the instantaneous ve- 
locity containing the Zitterbewegung, see Refs. [3, [Tg| . 
However, if one calculates the velocity using the en- 
ergy (34): Vi = de/dhki, it is the average velocity not 
depending on time and given by Eq. (25). This means 
that the two-band model in the matrix form still "keeps 
track" of the periodic Hamiltonian (1) from which it 
originates, because both give the Zitterbewegung. On 
the other hand, in the energy (34) the track of periodic- 
ity of the original Hamiltonian (1) is already lost. One 
can use the LK transformation to separate the conduc- 
tion and valence bands in the Hamiltonian (33), which 
gives e± = ±(e g /2 + ft 2 fc 2 /2toq) corresponding to the 
above mentioned expansion of the square root in Eq. (34) . 
In this case the velocity is vf — de±/dhki — ±hki/m,Q, 
i.e. it is the average velocity for each band with no Zitter- 
bewegung. Thus the two-band model is the simplest k ■ p 
description reproducing the essential features of the ini- 
tial periodic Hamiltonian (1). 

Finally, we want to consider briefly the important case 
of monolayer two-dimensional graphene in light of the 
above discussion. Graphene's band structure near the K 
point of the Brillouin zone is described by the Hamilto- 
nian 20] 



H = hi 



i . 



(41) 



where u w 10 8 cm/s. The above form can be considered 
to be a special case of the two-band model (33) with 
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the vanishing gap e g = and properly chosen matrix ele- 
ments irf 2 an d 7rf 2 - The resulting energy dispersion is lin- 
ear in quasi-momentum: e = ±uhk, where k = + k 2 . 

In view of our semi-relativistic analogy this case can be 
considered to be the "extreme relativistic limit". The 
matrix velocity operator Vi = dH/dhki does not com- 
mute with the Hamiltonian (41) and the instantaneous 
velocity contains the ZB component The veloc- 

ity Vi = de/dhki = uki/k represents an average velocity 
calculated in Eq. (21). The absolute value of velocity vec- 
tor for any direction is V = y ^ v 2 + v 2 = u. The velocity 

mass can still be defined as before: 1/m* = (l/h 2 k)de/dk. 
For the linear band dispersion one has de/dk — u, so 
that to* = Hk/u — e/u 2 . This gives, as before, e = m*u 2 , 
see Eq. (38). One can also write to* = e/u 2 which means 
that at the band edge (called in the literature " the Dirac 
point") the effective mass is zero, but as the energy in- 
creases the mass increases as well. Now let us suppose 
that an external force is applied along the x direction 
to an electron characterized by k y — 0. According to 
Eq. (12) there is d{hk x )/dt = F£ x . Since for k y = there 
is hk x = m*v x = m*u. Thus the change of hk x due to the 
external force goes entirely into the change of the mass. 



IV. VELOCITY AND ACCELERATION 
EFFECTIVE MASSES 

In this section we consider the use of the velocity ef- 
fective mass in spherical and spheroidal energy bands. 
At the end we mention some properties of the accelera- 
tion effective mass of charge carriers. We begin with the 
velocity mass to* which, as mentioned above, is much 
more useful than the acceleration mass. In our consider- 
ations below we are concerned with the average electron 
motion related to the quasi-momentum, so we drop the 
sign of "average" over the velocity, i.e. we write v = v. 
The important property of the velocity mass is that it 
is measured in the cyclotron resonance (CR). We first 
demonstrate it using the classical electron motion in a 
magnetic field. The equation of motion is 



d(hk) 
dt 



e(v x B), 



(42) 



where B = [0, 0, B] is a magnetic field applied along the z 
direction. Using the definition (20) of the velocity mass 
for a spherical band: hk = m*v, we have 



rn 



. dv 
~dt 



e(v x B). 



(43) 



Since, as we showed above, to* for a nonparabolic band 
depends in general on electron energy, one can imagine 
that m* depends also on time if the energy during the 
motion as not constant. However, it is well known that 
a constant and uniform magnetic field does not do any 
work, so the electron energy is a constant. For this rea- 
son we assumed rn* not to depend on time in arriving at 



Eq. (43). For the first two components the above equa- 
tion gives 



dvx 
~dt 

dVy 

~dt 



eB 
1 

TO* 

eB 



(44) 
(45) 



One can now differentiate Eq. (45) with respect to time, 
insert the result into Eq. (44) and arrive at the second- 
order differential equation for v v , which can be easily 
solved in terms of trigonometric functions. Instead, we 
simply guess the solutions (see Ref. [ll|): v x = vo cos(wf) 
and v y = — fosin(wt), in which to is the cyclotron fre- 
quency with which the electron circles on the orbit. Using 
the above solutions one obtains from Eqs. (44) and (45) 
the same result 



UJ 



eB 

TO* 



(46) 



Thus the cyclotron frequency is determined by the ve- 
locity mass to*. The cyclotron orbit can be obtained 
by integrating the velocity over time, which gives: (x — 
x o) 2 + (y — Vo) 2 — Vq/u 2 . This means that the classical 
cyclotron radius is also determined by to* . The same re- 
sult concerning the velocity mass can be obtained from 
the quantization of motion in a magnetic field. To be spe- 
cific, we use the orbital and spin quantization resulting 
from the band structure of InSb-type III-V semiconduct- 
ing compounds. The band structure includes three levels 
at the r point of the Brillouin zone (eight bands includ- 
ing spin). The resulting quantized orbital and spin levels 
are EHE?! 



e(n, k z , ±) = 



1/2 



where 



fiu>° 



± -q*u B B. 

2mS 2 yoP 



(47) 



(48) 



Here w° = eB/m*,, in which m*> is the band-edge mass 
[see Eq. (34)], <?q is the band-edge spin Lande factor 
and /is is the Bohr magneton. The cyclotron energy hu 
is given by the energy difference between two consecutive 
orbital levels: e(n+l,fc z ,±) — e(n,fc Z) ±). Using Eqs. (47) 
and (48) we have 



e(n + l) 2 - e(nf 
e(n + 1) + e(n) 



ht g eB 



1 



TOq e(n + 1) + e(n) ' 



(49) 



For small magnetic fields there is e(n + 1) + e(n) pa 2e, 
so that Huj « heB / '(2toqc/ 'e g ) — HeB/m*, see Eq. (35). 
Thus, again, the cyclotron frequency is determined by 
the velocity effective mass to*. The same reasoning can 
be applied to the so called "inverted " band structure 
of zero-gap and narrow-gap II- VI compounds based on 
HgTe and HgSe, see [H. 
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As mentioned above, the basic relation for the classical 
transport theory is hk = m*v. This leads to the defini- 
tion of carrier's mobility /i = qr/m*, which involves the 
relaxation time r and the velocity mass m* [24=] . As a 
consequence, the mobility is directly affected by the en- 
ergy variation of m*(e). Some d.c. transport phenomena 
at high magnetic fields do not depend on the relaxation 
time so that, by studying them, one gains a direct ac- 
cess to the mass m* [24|. Finally, the free-carrier optics 
depends on the band structure only through the veloc- 
ity mass. And so the reflectivity depends on (1/m*), 
the magneto-reflectivity on (l/m* 2 )/(l/m*), the Fara- 
day rotation is proportional to {1/m* 2 ) and the Voigt 
phase shift to (1/m* 3 ). Here the brackets denote appro- 
priate averages over electron energies in the band [24j . 
The knowledge of m* (e) gives then a direct information 
on the band structure. 

Another strong indication, that the velocity effective 
mass is much more useful than the acceleration mass, 
is the fact that the corresponding mass is commonly 
used in the special theory of relativity (STR). In STR 
this mass is defined by the relation: p = m{v)v, it is a 
scalar and it has the famous velocity dependence: m(v) = 
mo/(l — v 2 /c 2 ) 1//2 . Its energy dependence is not written 
down so often, but it is not difficult to derive. Since in 
STR the energy is given by E — [(m c 2 ) 2 +p 2 c 2 } 1 / 2 , the 
velocity is Vi = dE/dpi — piC 2 /E, and the velocity mass 
is m = E/c 2 = 2moE / (2ni()C 2 ) . It is seen that, using 
the semi-relativistic analogy: mo — > m*, and 2moc 2 — > e g , 
the relativistic velocity mass has the same energy depen- 
dence as the effective velocity mass resulting from the 
two- band k ■ p model, Eq. (35). We mention that the 
relativistic velocity-dependent mass m(v) is somewhat 
reluctantly used in STR by some authors because of its 
unorthodox transformation properties, see 26]. However, 
the transformation problem is not relevant for solids. 

To conclude our considerations of the velocity mass 
we treat an important case of ellipsoidal energy bands 
which occur in semiconducting II- VI lead salts PbTe, 
PbSe, PbS, as well as in silicon and germanium. Such 
an energy band with arbitrary nonparabolicity can be 
described by the relation [U, [2{| 

7( e ) = aapkakp, (50) 

where 7(e) is a "reasonable" function of energy describ- 
ing the nonparabolicity of the band. The limiting as- 
sumption is that the shape of the ellipsoid does not vary 
with the energy. We use the sum convention over the re- 
peated coordinate indices. The tensor a a /3 is symmetric 
and it can be brought to a diagonal form by an appro- 
priate rotation of coordinates in k space. Then the un- 
equal diagonal components an express band's ellipsoidal 
shape. An inverse tensor of velocity mass is defined by 
the relation (22). On the other hand, since the velocity 
is Vi = de/dhki, one obtains with the help of Eq. (50) 



For a parabolic band there is 7(e) = e, and the in- 
verse mass components are numbers. A spherical band 
with arbitrary nonparabolicity is described by a a /3 = 
S a p, the Kronecker delta. This gives 7(e) = k 2 , [see 
Eq. (50)], which is equivalent to the spherical case con- 
sidered above. 

Finally, we calculate the inverse tensor of acceleration 
mass for a spherical nonparabolic band e(fc). The general 
expression of the inverse mass tensor relating force to 
acceleration is well known 

(J_\ 1 d 2 e 

n 2 dkidkj- 1 ' 

A simple manipulation gives 

(J_\ = _L^X h^fj 2 ^ _L^\ (ro, 

Uv« h 2 kdk v k 2 \h 2 dk 2 h 2 kdk) y ' 

It is seen that, in contrast to the velocity mass, the accel- 
eration mass is not a scalar quantity even for a spherical 
energy band. Interestingly, it is the band nonparabolic- 
ity that makes the acceleration mass non-scalar. For a 
standard parabolic band: e = H 2 k 2 /2rriQ, the second term 
in Eq. (53) vanishes, so that (1/M*)y = (l/m%)5ij. This 
is identical with the velocity mass for a parabolic band. 
Only in this simple case the two masses coincide. Because 
of the non-scalar character of the acceleration mass (53) 
the acceleration in a nonparabolic band is not parallel to 
the force. This feature is well known in the special rel- 
ativity, which illustrates once again the semi-relativistic 
analogy. 



V. CONCLUSIONS AND SUMMARY 

We summarize our work by enumerating the main con- 
clusions and indicating the corresponding equations. For 
a carrier moving in a periodic potential, the momentum, 
velocity, and acceleration are not constants of the mo- 
tion, see Eqs. (l)-(3). The quasi-momentum, which is 
a distinctly different operator from the momentum, see 
Eq. (8), is a constant of the motion in a Bloch state, see 
Eq. (9). The average electron velocity in a Bloch state 
is given by a gradient of the energy with respect to the 
quasi-momentum hk, see Eq. (19). It is this average ve- 
locity v which is used in the classical transport theory 
for charge carriers. A "velocity effective mass" is defined 
as a quantity relating the average velocity to the quasi- 
momentum, see Eqs. (20) and (21). The velocity mass 
for a spherical energy band is a scalar, see Eq. (23), and 
it enters into the basic relation for the transport theory, 
see Eq. (25). A two-band k ■ p model in the matrix form 
is the simplest description of the band structure that still 
keeps track of the periodic potential, see Eq. (33). The 
two-band k ■ p Hamiltonian (33) and the resulting en- 
ergy (34) bear strong similarity to the description of free 
relativistic electrons in a vacuum. In particular, they 
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lead to an analog of the famous Einstein relation be- 
tween the mass and the energy, see Eq. (38). In this 
perspective, the band structure of gapless graphene can 
be regarded as an extreme rclativistic case. The velocity 
effective mass is much more useful than the acceleration 
mass commonly introduced in solid state textbooks. In 
particular, it is the velocity mass that is measured in 
the cyclotron resonance, see Eqs. (46) and (49), in d.c. 
transport phenomena and in the free-carrier optics. The 



velocity mass can also be introduced for ellipsoidal non- 
parabolic energy bands, see Eq. (51). 
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